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Electronic structure and dielectric property in an electronic ferroelectricty, where electric polarization 
is driven by an electronic charge order without inversion symmetry, are studied. Motivated from layered 
iron oxides, roles of quantum fluctuation on ferroelectricity in a paired-triangular lattice are focused on. 
Three types of the extended V — t model are examined by the variational Monte-Carlo method with the 
Gutzwiller-type correlation factor. It is shown that electron transfer between the triangular layers corre- 
sponding to the inter-layer polarization fluctuation promotes the three-fold charge order associated with 
an electric polarization. This result is in highly contrast to the usual manner of quantum fluctuation in the 
hydrogen-bond type ferroelectricities and the quantum paraelectric oxides. Spin degree of freedom of elec- 
tron and a realistic crystal structure for the layered iron oxides further stabilize the polar charge ordered 
state. Implications of the numerical results for layered iron oxides are discussed. 

KEYWORDS: Ferroelectricity, Frustration, Charge Order 



1. Introduction 

Roles of electronic degree of freedom in ferroelectric- 
ity and related phenomena have attracted much attention 
for a long time. In particular, co valency contribution on the 
displasive-type ferroelectricity has been revealed by the mod- 
ern afc-initio electronic-structure calculation and the Berry 
phase theory for the electric polarization. Recently discov- 
ered multiferroics, i.e. coexistence of ferroelectricty and mag- 
netism, are another example. 1_6) Through the intensive theo- 
retical and experimental examinations, it has been uncovered 
that the symmetric/anti-symmetric exchange interactions and 
the exchange striction effect are origin of this type ferroelec- 
tricity. 

It is known that there is different type of ferroelec- 
tricity driven by electronic degree of freedom. A long- 
range order of the electronic charge without inversion sym- 
metry is responsible for a macroscopic electric polariza- 
tion. This phenomenon is termed the electronic ferroelec- 
tricity or the charge-order type ferroelectricity. 7 ' 8) A pos- 
sibility of this type of ferroelectricity has been suggested 
experimentally in a variety of transition-metal oxides, e.g. 
Pr(Sr a: Cai_ a; )Mn 2 07, 9) LuFe 2 4 , 10 ' n > low-dimensional or- 
ganic salts, e.g. TMTTF 2 X (X=PF 6 , AsF 6 ), 12) a-(BEDT- 
TTF) 2 I 3 , 13) k-(BEDT-TTF) 2 Cu 2 (CN) 3 14 - 16) and others. Most 
of this type of material belongs to the so-called quarter-filled 
system where the number of carrier per site is 0.5. The elec- 
tronic ferroelectricity has some similarities to the hydrogen- 
bond type ferroelectricity; localization and derealization of 
charged particles, i.e. electron or proton, are concerned in the 
ferroelectric transition. One noticeable difference between the 
two is that the electron mass is much smaller than the pro- 
ton mass. Therefore, large dielectric fluctuation is expected to 
play crucial roles in the electronic ferroelectricity. 

Layered iron oxide i?Fe 2 04 (i?=Lu, Yb, Y) is recog- 
nized as a representative electronic ferroelectricity. The crys- 
tal structure in i?Fe 2 04 consists of a stacking of the R-O 



layers and the Fe-O layers where Fe ions form the paired- 
triangular lattices. 17 ' A nominal valence of an Fe ion is 2.5+, 
and an equal amount of Fe 2+ with a d 6 (S = 2) configuration 
and Fe 3+ with d (S = 5/2) coexists. A charge order of Fe 
3d electrons in LuFe 2 04 was found by the x-ray and electron 
diffraction experiments 1 8-2 where the superlattice peaks ap- 
pear at (n/3 n/3 3m + 1/2) below about 320K. A ferrimag- 
netic long-range order was confirmed by the neutron diffrac- 
tion experiments where the magnetic peaks at (1/3 1/3 m) 
are observed below about 250K. 22 ~ 24) The temperature de- 
pendence of the electric polarization was measured by the 
pyroelectric current. It was found that the polarization ap- 
pears around 320K and increases around 250K. 10 ' n) These 
results are interpreted that the electric polarization is deduced 
by the Fe 2+ /Fe 3+ charge order, and is strongly coupled with 
the spins in Fe ions. 

A primitive model for the observed charge order in i?Fe 2 04 
associated with the electric polarization was first proposed 
by Yamada and coworkers. 18, 19) It was considered that Fe 
ions are aligned as • • Fe 3+ -Fe 3+ - Fe 2+ • • • along one of 
the crystal axes in a upper triangular layer and as ■ • • Fe 3+ - 
Fe 2+ -Fe 2+ • ■ • in a lower layer. As a result, charge imbalance 
in a paired triangular lattice induces the electric dipole mo- 
ment. Recently, detailed microscopic calculations for stabil- 
ity of the polar CO structure have been presented. 25 ~ 30) One of 
the present authors and coworkers presented a microscpic the- 
ory for the electronic structure and dielectric properties. 26 ' 27 -* 
They treated the electronic charge and spin degrees of free- 
dom as classical variables in the calculation, and showed that 
the polar charge order does not appear at zero temperature 
and is stabilized by the thermal fluctuation effect. This result 
is consistent with the recent optical experiments where large 
charge fluctuation remains far below the charge ordering tem- 
perature. 31 ' 

As well as the thermal fluctuation, the quantum fluctuation 
of the electronic charge, i.e. the electron transfer between the 
sites where the electronic potential takes its minima, is ex- 
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pected to have major contribution to this-type of ferroelectric 
transition. This is because most of this type of ferroelectric - 
ity is located at vicinity of the metal-insulator transition. Role 
of the quantum fluctuation on ferroelectricity has been exam- 
ined for a long time in the hydrogen-bond type ferroelectric- 
ity, 32 ' 33 ' i.e. the proton tunneling between the double-potential 
well, and the quantum paraelectric oxides, such as SrTiOs 
and KTaC>3. 34) It is known that the quantum motion of pro- 
tons/ions prevents the system from the long-range ferroelec- 
tric order. More remarkable quantum fluctuation effect is ex- 
pected in the electronic ferroelectricity, because the electron 
mass is much smaller than the ion/proton mass. It is highly 
nontrivial whether the conventional quantum fluctuation ef- 
fects seen in the hydrogen-bond type ferroelectricy and the 
quantum paralectric oxides are naively applicable to the elec- 
tronic ferroelectricity or not. 

The present paper addresses the issue of the quantum fluc- 
tuation effect in the electronic ferroelectricity, in particular, 
the dielectric and magneto-dielectric properties in the lay- 
ered iron oxides. We introduce the three-types of the extended 

V — t model in a paired-triangular lattice and focus on the 
electron transfer effect between the layers, corresponding to 
the quantum fluctuation of the electric dipole moment. The 
models are analyzed by the variational Monte-Carlo (VMC) 
method by which the quantum fluctuation effect, the long- 
range Coulomb interaction and the geometrical frustration are 
treated properly. It is shown that the electron transfer between 
the layers tends to stabilize the polar charge order, in con- 
trast to the conventional manner of the quantum fluctuation 
effect on the ferroelectricity. The spin degree of freedom and 
the realistic lattice structure in i?Fe2C<4 reinforce the electric 
polarization. 

In Sect. II, the extended V—t model in the paired-triangular 
lattice and the formulation of the VMC method are intro- 
duced, and the numerical results in this model are presented. 
In Sect. Ill, roles of the spin degree of freedom in the di- 
electric and magneto-dielectric properties are examined. In 
Sect. IV, the numerical results in the model, where the real- 
istic crystal structure in i?Fe2C>4 is taken into account, are 
presented. Section V is devoted to summary and discussion. 
A part of the present results were briefly reported in Ref. 35) 

2. Electric Polarization in V — t Model 

As mentioned in the previous section, the crystal lattice of 
i?Fe2C<4 consists of the alternate stacking of the Fe-O double- 
triangle layers, termed W-layer, and the R-O layers along the 
c axis. A unit cell with the R3c symmetry includes the three 
W-layers and the three R-O layers. We focus on electronic 
structure in a W-layer which dominates the electric and mag- 
netic properties in i?Fe204. We consider the three types of 
the V — t model in the W-layer, and start from the spin-less 

V — t model in a regularly stacked paired-triangular lattice. 
The double triangular lattices are stacked along the c axis as 
shown in Fig. 1 . The long-range Coulomb interactions and the 
inter-site electron transfers between the spin-less fermions are 
taken into account. A way of a stacking of the triangular lat- 
tices is somewhat different from that in i?Fe2C>4. Effects of 
the spin degrees of freedom and the realistic stacking of the 
triangular layers will be introduced in Sects. 3 and 4, respec- 
tively. 
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Fig. 1 . (Color online) A pair of the triangular lattices and interactions in 

H V t. 



The model Hamiltonian is given as 
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where & ■ is the creation operator for a spin-less fermion at 
site i on the upper layer (m = u) or the lower layer (m = I), 
c m i is the number operator. A subscript i takes 



and n, 



a two-dimensional coordinate in a triangular lattice. To exam- 
ine roles of frustration, we consider an anisotropic triangular 
lattice where bonds along the [110] direction are inequivalent 
to the bonds along [100] and [010] as shown in Fig. 1. We 
introduce the three kinds of the transfer integrals t, t' and t± 
for the [100]/[010], [110] and [001] directions, respectively. 
In the same way, we introduce the three kinds of the inter-site 
Coulomb interactions V, V and V± for the nearest-neighbor 
(NN) pairs. Symbols an< ^ represent summations 

for the NN pairs in the [100]/[010] and [110] directions, re- 
spectively. Ratios V jV and t! /t represent magnitude of frus- 
tration. The number of the fermion per site is 0.5. 

We analyze this model by using the VMC method. A vari- 
ational wave function is given by a following product form, 



(4) 



where $ is the one-body Hartree-Fock (HF) part and V is the 
many-body correlation factor. We introduce the Gutzwiller- 
type correlation factor given by 

V = (l-vn mi n mj ) ] [ (l-v'n mi n m j) 

m(ij) 



m<ij> 



1 (1 - v ± n ui nii), 



(5) 



where v, v' and v± are the variational parameters. As shown 
in Fig. 2, in the one-body part of the wave function, we con- 
sider the following four-types CO structures: (i) a two-fold 
CO along the [110] direction termed COi/- 2 -l, (ii) a two-fold 
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(a) CO 




(b) CO 



(c) CO 




Fig. 2. (Color online) Schematic CO structures cosidered in Wvt- (a) 
C0 1/2 -I, (b) C0 1/2 -II, (c) CO1/3 and (d) C0 1/4 . The inset of (c) is a 
side view along the [120] direction. 



CO along [100] termed COi/ 2 -II, (iii) a three-fold CO along 
[110] termed CO1/3, and (iv) a four-fold CO along [110] 
termed CO1/4. These CO structures are the candidates for 
the mean-field (MF) ground state. In the COx/2-F CO^-II 
and COi/4 structures, stackings of the CO patterns in the two 
triangular-layers are out-of-phase. In each CO, the one-body 
part of the wave function is given by the HF wave function 
where the corresponding CO's are assumed in the MF order 
parameter. We take that black, white and grey circles in Fig. 2 
correspond to the MF charge densities of (n TO , ) — 1/2 = A, 
—A and zero, respectively. The MF order parameter A is re- 
garded as a variational parameter and is optimized in $. For 
example, $ for the CO1/3 structure is obtained by diagonal- 
izing the following HF Hamiltonian, 
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v-t — 
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We define 

T k = -t{e lkl + e lk2 ) - t' e - i[kl+k2 \ 

and 

W = -A(2V + V + V±_), 

where k\ and k-2 are the wave-vectors along the [100] and 
[010] directions, respectively. A subscript A(= A, B, C) rep- 
resents the three sublattices in the CO^ structure shown in 
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Fig. 3. (Color online) Phase diagram in Hvt< The relations t = t' = tj 
and V = V± are imposed. 



Fig. 2(c), and the summation ^ k runs over the reduced two- 
dimensional Brillouin zone. The operator c m A k is the Fourier 
transform of CmXi where site coordinate is redefined by the 
two subscripts A and i. The one-body parts of the wave- 
functions for other CO's are defined in the same way. A po- 
lar CO1/3 structure of the present interest has a possibility to 
be realized in the CO1/3 structure, when all the grey circles 
in the upper (lower) plane in Fig. 2(c) becomes charge rich 
(poor) sites due to the correlation factor V in a case of the full 
electric polarization. 

In the numerical calculation, we adopt 2 x 10 5 — 5x1 5 
samples in most of the VMC simulations. The fixed-sampling 
method is used to optimize the variational parameters. 36) An 
error of the energy expectation value is of the order of 10 -4 £. 
Cluster sizes are of N s = L x L x 2(= 2JV) (L < 12) sites 
with the periodic-boundary condition and the antiperiodic- 
boundary condition. 

In Fig. 3, the phase diagram is presented in the plane of 
V'/V and t/V. We chose t = t' = t x and V = V x . This 
phase diagram is obtained from the calculated energy and 
charge correlation function defined by 

1 



A/(q) 



ijm 



,iq-(R„ 



-R„ 



,)-n 2 ), (12) 



where n = 1/2 and R m i is a position of site i on the m layer. 
In the classical limit, i.e. t = 0, the CO1/2-I and CO1/4 struc- 
tures are realized in the regions of V'/V < 1 and V'/V > 1, 
respectively. At V'/V = 1 with t = 0, three CO structures, 
i.e. CO1/2-II, CO1/4, and CO1/3, are degenerate. When the 
electron transfer t is introduced, the COi/3 structure is stabi- 
lized in a finite parameter region. 

Let us focus on the results in V'/V = 1 in more detail. 
In Fig. 4(a), the energy expectation E = N^ 1 (Hvt) is plot- 
ted as a function of t/V. Reduction of E in CO1/3 with in- 
creasing t/V is more remarkable than those in the other CO 
structures. A stabilization of the COiy 3 structure is attributed 
to the ^-linear dependence of E, in contrast to the quadratic 
dependence in other CO's. The charge correlation function 
at q = (27r/3, 2tt/3)(= qi/3) is presented in Fig. 4(b). A 
possible maximum value of N(qx/3) is 1/97V for the CO1/3 
structure. Magnitude of N(qi/ 3 ) is about 70% of its maxi- 
mum value below t/V = 0.4 and is almost independent of 
the system size L. This correlation almost disappears above 
t/V = 0.4 where the metallic phase appears as shown later. 
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Fig. 4. (Color online) (a) Energies for several CO's. Solid and broken lines 
are for the energies obtained by the second-order perturbation in CO1/3 
and other CO's, respectively, (b) Charge correlation function iV(q) at q = 
(27r/3, 2n/3). (c) Polarization correlation P. Relations t = t' = tj_ and 
V = V = V± are imposed. Open and filled symbols are for the results 
in L = 6 and 12, respectively. 
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Fig. 5. (Color online) Polarization correlation P. Parameters are chosed 
to be V'/V = 0.95 in (a) and 1.05 in (b). Relations t = t' = t ± and 
V = V± are imposed. Open and filled symbols are for the results in L = 6 
and 12, respectively. 



In Fig. 4(c), we plot the correlation function of the electric 
polarization defined by 



P=(p 2 )\ 



(13) 



(14) 



A possible maximum value of P in CO1/3 is one. We also 
present the results for V'/V = 0.95 and 1.05 in Fig. 5 to 
compare P's in several CO's. Magnitude of P in COx/3 is 
much larger than those in CO^-I and COi/$. However, val- 
ues of P decrease almost by half with increasing the system 
size L from 6 to 12. These results indicate that the electric po- 
larization in the thermodynamic limit is supposed to be much 
smaller than the maximum value. We conclude that a large 
polarization fluctuation appears in the CO1/3 phase, in com- 
parison with other CO's, but the robust electric polarization is 
not expected in the thermodynamic limit. 

As explained below, the energy gain in COx/3 is caused by 
the inter-layer electron transfer, t±. Let us focus on a pair of 
sites represented by grey circles in the upper and lower planes 
in Fig. 2(c). Since these sites are surrounded by the in-plane 
three NN charge-rich sites and the NN three poor sites, the 
Coulomb interactions between NN sites are canceled out, and 
the classical energy does not depend on the charge configu- 
rations at the pair of grey-circle sites. Thus, the stabilization 
of the COi/3 structure is caused by the first order of t± be- 
tween the grey-circle sites. This situation does not realize in 
other CO structures. Even when the longer-range Coulomb 
interactions than V, V and Vj_ are taken into account, the 
above scenario for stabilization of CO1/3 is expected to sur- 
vive in the case where t± is larger than those interactions. 
To support this scenario, the energy expectations are calcu- 
lated in the second-order perturbation with respect to t for 
several CO's. The results are shown in Fig. 4(a). In the region 
of small t, the i-linear dependence of E in COx/3 obtained 
by the VMC method is well reproduced by the perturbational 
calculation. It is concluded that the COiy3 is stabilized by a 
combined effect of the the geometrical frustration effect in the 
triangular lattice and the electron transfer between the layers. 

We briefly touch the metal-insulator transition in this 
model. Discontinuous changes in P and A r (q 1 / 3 ) at t/V ~ 
0.4 shown in Figs. 4(b) and 4(c) are associated with the metal- 
insulator transition. This is directly confirmed by calculating 
the momentum-distribution function defined by 



n(k) 



_L ik (Ri-Rj 



(15) 



mi j 



In Fig. 6, n(k) at V'/V = 1 is shown for various values of 
t/V along the high-symmetry lines in the Brillouin zone. In 
t/V > 0.4, discontinuous changes are seen in n(k) along 
(0,0)-(tt,0) and (tt, tt)-(0, 0). 

Let us compare the present results with the ones in the 
single-layer V — t model. It was reported in Refs. 37 ~ 39) that 
a long-range three-fold CO coexists with a metallic character 
in the case of V = V and t = t' in the single-layer V — t 
model. This is termed a "pinball liquid" (PBL) state. Here we 
set up the trial wave function for the PBL state in the paired- 
triangular lattice obtained by diagonalizing the following HF 
Hamiltonian, 



p 1/3 
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Fig. 6. (Color online) The momentum distribution function n(k) for var- 
ious values of t/V. Relations t = t' = t± and V = V = V± are 
imposed. System size is L = 12. 




0.4 0.6 

VJ V 

Fig. 7. (Color online) Energies in the hole-pinned pinball liquid state and 
the COj/3 state. Relations t = t' = tx and V = V are imposed. Open 
and filled symbols are for the results in L = 6 and 12, respectively. 
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We introduce ki and hi as the wave-vectors along the [100] 
and [010] directions, respectively, and define 

1, 



Wf = - 7p A(2V + V + -V±), 



and 



v±), 



V-ivj, 



(19) 



(20) 



(21) 



2 2 

where j p is a numerical factor taking 1 and — 1. When j p = 
1, the trial wave function represents the so-called electron- 
pinned PBL state where one of the three sublattices in the 
triangular lattice is occupied by 1/3-electrons per site and re- 
maining 1/6-electrons per site move on other two subllatices. 
In this state, the mean fields in three sublattices are taken to 
be (n ui ) - 1/2 = (A, -A/2, -A/2) and (n H ) - 1/2 = 
(— A/2, A, —A/2). When 7 p = —1, the wave function rep- 
resents the so-called hole-pinned PBL state, where one of the 
three sublattices is occupied by 1/3-holes per site and 1/6- 
holes per site move on others. The mean fields are taken 
to be (n ui ) - 1/2 = (-A, A/2, A/2) and (n H ) - 1/2 = 
(A/2, -A, A/2). 

The stabilities for the two PBL states are compared with 
the COx/3 state. The electron-pinned PBL has much higher 
energy in comparison with COx/3 in the present parameter 
range. In Fig. 7, we present the energy expectations of the 
hole-pinned PBL state and the CO!/ 3 state att = t' = t± and 
V = V . The PBL is realized only in a region of the small 
inter-layer Coulomb interaction, V±/V ^ 0.1, where the two 
triangular layer is almost independent. This is consistent with 



the previous results in the single triangular layer V — t model 
analyzed by the VMC method. 39 ' 

3. Role of Spin Order 

So far, we examine the spin-less V — t model in a paired 
triangular lattice. It is shown that the inter-plane transfer in- 
tegral stabilizes the COi/ 3 phase where a larger polarization 
fluctuation appears in comparison with other types of CO's. 
However, magnitude of the electric polarization in the ther- 
modynamic limit is supposed to be much smaller than the 
possible maximum value. We consider that additional small 
factors, which are not taken into account so far, may stabi- 
lize the electric polarization in the COx/3 phase. In this sec- 
tion, we investigate roles of the spin degree of freedom on 
the dielectric and magneto-dielectric properties in a paired- 
triangular lattice. 

We consider the spin degree of freedom on the basis of 
LuFe2C»4 where two kinds of spin, 5 = 2 for Fe 2+ and 
5 = 5/2 for Fe 3+ , exist. The electron configurations of Fe 2+ 
and Fe 3+ are 3d 6 and 3d 5 , respectively. It is thought that an 
excess electron in Fe 2+ occupy one of the doubly degener- 
ate d x 2_ y 2 and d xy orbitals which have a larger transfer inte- 
gral in the plane than other 3d orbitals. 26 27) We introduce, in 
the model Hamiltonian, a localized spin with 5 = 2 at each 
site, itinerant carriers corresponding to the excess electrons 
in Fe 2+ , and the ferromagnetic Hund couple between them at 
each site. The number of itinerant carriers per site is set to be 
0.5. This is the generalized double-exchange model. For sim- 
plicity, we assume that the Hund coupling is infinite, and the 
localized spins are the Ising type denoted by a m i (m = u, I) 
which takes ±1. The Hamiltonian is given by 



Hafm- 



(22) 



The first and second terms represent the electron transfer and 
the inter-site Coulomb interaction, respectively, given by 
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<2j>ms 
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Fig. 8. (Color online) Interactions in %g. Large and small arrows repre- 
sent spin directions of the localized spins and the conduction electrons, 
respectively. 
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(24) 



where c m i S is the electron annihilation operator at site i on the 
m layer with spin s(=t: i)- The third term in Eq. (22) is for 
the antiferromagnetic superexchange interaction between the 
localized spins given by 



H A 



I'M 



<ij>m f» ™ 

E 



(ij)m 



(25) 



where the bond-dependence of the antiferromagnetic superex- 
change interaction J$ is neglected. In Hf , we consider that 
the transfer integrals depend on configuration of the localized 
spins. In the infinite-limit of the Hund coupling, it is reason- 
able to assume the following form of the transfer integral, 



to 



t, 
0, 



for <t„ 
for o ri 



(26) 



7^ ^mj j 

where t is introduced in Eq. (2). We also define ^ mj(Tm . and 

t±a m iCTm- i R ?H i R *h e same wa Y- A schematic view of the 
Hamiltonian is presented in Fig. 8. As for the localized spin 
structure, we consider all possible ordered structures char- 
acterized by the momenta (2tt/3,2tt/3) which corresponds 
to the observed neutron diffraction peak at (1/3 1/3 m) in 
LuFe2C>4. 22 ~ 24) These magnetic structures are shown in Fig. 9 
where equivalent patterns are not shown. 

We analyze this model Hamiltonian by using the VMC 
method. The one-body part of the variational wave function 
$ introduced in Eq. (4) for the CO1/3 structure is obtained 
from the HF Hamiltonian. Here, we redefine the HF Hamilto- 
nian in the three sublattices as follows, 

h s 1/3 = E 

km s 

-t± E ), (27) 

kAs 



with 








(28) 



"km 



la mA (T m 





where we introduce ~f a 



= S n 





and 77 auA(Ti 
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The correlation factor V is defined in Eq. (5). The en- 









Fig. 9. (Color online) Configurations of the localized spin direction con- 
sidered in H s . 
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ergy is optimized in each spin configuration shown in Fig. 9. 



Fig. 10. (Color online) Energy optimized in each spin configuration in H s . 
Numbers ( l)-(6) correspond to the spin configurations given in Fig. 9. Pa- 
rameters are chosen to be t/V = 0.2, t = t' = tj_ and V — V = V±. 



A ratio of the number of the spin-up and spin-down electrons, 
N+/Nu is also optimized. We restrict our calculations to the 
case of V = V = V ± and t = f = t±. 

In Fig. 10, the optimized energies in the six spin config- 
urations are plotted as a function of J$- It is shown that the 
energy in the spin configuration 4 (see Fig. 9) is the lowest. In 
Fig. 11, the energy and the electric polarization correlation in 
the spin configuration (4) are plotted as a function of N^/N e 
where N e = 7V t + N±. At N t /N e = 0.33, the energy minima 
and a large magnitude of the polarization are realized. These 
magnitudes are almost independence of the system size L. 
That is, the spin degree of freedom induces the robust electric 
^polarization. 

The spin structure realized at the energy minima shown in 
Fig. 1 1(a) is presented in Fig. 12. Let us focus on sites A and 
B in the upper layer and consider the electron transfer from 
these sites to the sites just below them. When the spin degree 
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Fig. 1 1 . (Color online) (a) Energy and (b) polarization correlation as func- 
tions of the number of the up-spin conduction electrons in H s . Param- 
eters are chosen to be t/V = 0.2, Js/t = 0.05, t = t' = tj_ and 
V = V = V± . Open and filled symbols are for the results in L = 6 and 
12, respectively. 



[120] [110] 
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Fig. 13. (Color online) Lattice structure and interactions in Ww- In 
i?Fe204, positions of the oxygen ions in the lower-layer are indicated by 
small black circles. 



4. Effect of Realistic Crystal Structure in HFe 2 4 

In this section, we introduce the realistic crystal structure 
in i?Fe2C>4 as a candidate to stabilize the electric polariza- 
tion in the CO1/3 structure. In i?Fe2C>4, a position of a Fe 
ion in the upper layer is not upon a Fe ion in the lower layer, 
but is upon an O ion surrounded by three in-plane Fe ions as 
shown in Fig. 13. We modify the Hamiltonian Hvt defined 
in Eq. (1) by taking this crystal structure into account. The 
model Hamiltonian is given by 




BR 



IT 



ftA [10 0] 



Fig. 12. (Color online) Left: spin structure in the optimized state corre- 
sponding to the energy minimum at Nf/N e = 0.33 in Fig. 11(a). Large 
and small arrows represent spin directions of the localized spin and the 
conduction electrons, respectively. Right: a side view along the [120] di- 
rection. 



of freedom is neglected, the classical energy is not changed by 
these electron transfer processes. This is the similar situation 
to the COi/3 structure in the spin-less V — t model as ex- 
plained in the previous section. Therefore, the energy in this 
structure is gained by the linear order of the electron trans- 
fer between the layers. However, in contrast to the previous 
spin-less case, the spin order suppresses the electron transfer 
between the layers and promotes the charge imbalance be- 
tween the two layers as follows; 1) In the upper layer, 2/3 of 
the localized spins are polarized toward the upper direction. 
Therefore, the up-spin conduction electrons (A) are confined 
in the upper layer to gain the in-plane kinetic energy. 2) The 
transfer of the down-spin electron (B) from the upper layer to 
the lower layer is prevented, because increasing of the num- 
ber of the down-spin electron in the lower layer enhances the 
kinetic energy. As the results, 2/3 of the electrons are confined 
in the upper layers and the charge imbalance between the lay- 
ers occurs. That is, the electric and magnetic polarizations are 
induced cooperatively. 



ft — ftf T fty, 



(30) 



with 



^ab c lni C ?nj t c c \ii c lj 

<ij>m (ij) 

+ H.C., (31) 
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A schematic view of the lattice structure and the interactions 
are shown in Fig. 13. We consider the long-range Coulomb 
interactions between the in-plane NN sites (V^nnX those be- 
tween the inter-plane NN sites (T4nnX and those between the 
inter-plane next-nearest-neighbor (NNN) sites (T4nnn)- Sym- 
bols Yl<ij>> m ^ Sfij] represent summations for the 
pairs corresponding to V^bNN, Knn and T4nnn, respectively. 
When the 1 / r-type Coulomb interaction is assumed, we have 
VcNN/KbNN = 1.2 and VcNNN/^abNN = 0.77 for the crystal 
structure in i?Fe204. 17) As for the electron transfer term "H^, 
the transfer integral between the in-plane NN sites (i a b), and 
that between the inter-plane NN sites (i c ) are considered. 

This model is analyzed by the VMC method. In the many- 
body correlation term V, we consider the variational param- 
eters t> a bNN, w C NN and Dcnnn which correspond to the in- 
teractions V'abNN, V^nn and T4nnn, respectively, instead of 
Eq. (5). In the one-body term of the wave function $, we 
adopt the following CO structures: (i) a three-fold CO along 
the [100] direction termed COiy 3 , (ii) a six-fold CO along 
the [100] direction termed COjyg, (iii) two types of two-fold 
CO's along [100], i.e. (iii) COi /2 -A and (iv) COi/ 2 -B, and (v) 
afour-foldCO along [100] termed CO1/4. These are schemat- 
ically shown in Fig. 14. We define that the CO1/3 structure 
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(a) polar CO 



(b) CO, 




Fig. 14. (Color online) Schematic CO structures considered in "Hy t '. (a) 
polar COi/3, (b) C0 1/6 , (c) C0 1/2 -A, (d) C0 1/2 -B and (e) C0 1/4 . The 
insets of (c) and (d) are side views along the [120] direction. 



in this realistic lattice has a polar CO configuration which is 
schematically given as---»o»»o»--- (•••o»oo»o---) along 
[100] in the upper (lower) plane. Symbols o and • are for the 
charge densities of (n m i — 1/2) = A and — A, respectively, 
where A is the variational parameter in the one-body part of 
the wave function. When this CO is fully ordered, the electric 
polarization p in Eq. (14) is one. The CO configuration along 
[100] in the non-polar CO^g structure is schematically given 
as ••••ottoo' - (•••o»oo»»-- - )in the upper (lower) 
plane. In each CO's, the one-body part of the wave function 
is obtained from the HF Hamiltonian. For example, the HF 
Hamiltonian for the polar CO1/3 is given by 



H R 1/3 
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Fig. 15. (Color online) Phase diagram in Hy t - Parameters are chosen to 
be t c /*ab = 1 and V cNN /V abNN = 1. 
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Fig. 16. (Color online) (a) Differences between the energies in several 
CO's and that in the CO]y3 obtained in Hy t . (b) Polarization correla- 
tion P. Parameters are chosen to be = t c and V a bNN = V c nn = 
2V c nnn- Open and filled symbols are for the results in L = 6 and 12, 
respectively. 
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The summation ^ k runs over the Brillouin zone for the CO 
structure with the momentum (2w/3, 2tt/3). The charge den- 
sities in sublattices A and B in the upper plane and that in C in 
the lower plane are A(> 0), and those in other sublattices are 
—A. The one-body parts of the wave functions in other CO's 
are also obtained in the same way. 

In Fig. 15, the phase diagram in the plane of T4nnn /VabNN 
and t a b/ VLbNN is presented. We chose the relations V^nn = 
VabNN and i a b = t c , although the results introduced below 
are robust qualitatively around these parameter. This phase 
diagram is obtained by the energy and the charge correlation 
function introduced in Eq. (12). A point of V C NNN/V a bNN = 
0.5 at tab/ VabNN = is a fully frustrated point which cor- 
responds to V jV = 1 and t/V = in Fig. 3. In this point, 
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(a) polar CO, 



(b) CO,, 





Fig. 17. (Color online) Inter-layer charge fluctuation processes by t c in the 
polar COx/3 phase (a) and those in the CO^g phase (b). The bold arrows 
represent the charge fluctuation where the energy reduction is of the order 
of tc/^cNNN- Large (small) circles represent the sites in the upper (lower) 
layer. 



the polar CO^ is degenerate with CO^ and COiy2-A. It 
is shown that the polar CO1/3 is realized in finite parameter 
region, when the electron transfer is taken into account. 

We show the results in V^nnn/Kann = 0.5 in more de- 
tail. In Fig. 16(a), we plot the energy differences of the several 
CO's from that in the COi/ 3 structure. The result for COi/ 2 - 
B is much larger than the scale in the figure. The polar COi/ 3 
phase competes mainly with the non-polar CO1/6, and is sta- 
bilized in fab/VabNN ^ 0.15. The electric polarization cor- 
relation is presented in Fig. 16(b). The results in the COiy 3 
phase are independent of the system size L, and show almost 
its maximum value. On the other hand, the value of P in the 
CO1/6 phase are much smaller than those in the CO1/3 phase 
and decrease with increasing of the system size. 

A stability of the polar CO1/3 is attributed to a combina- 
tion of the inter-layer charge transfer, t c , and the long-range 
Coulomb interaction, T4nnn, in this realistic crystal structure 
as follows. Let us focus on the bonds represented by the gray 
bars in the Fig. 14(a). These bonds connect the black-and- 
white circles and are responsible for the electric dipole mo- 
ments. The inter-layer Coulomb interaction Knnn causes the 
interaction along the [120] direction between the bonds, and 
induces the ferroelectric interaction between the dipole mo- 
ments along this direction. Although the dipole moments are 
induced in the local bonds, there is still competition between 
the polar CO1/3 and the antiferro-electric CO1/6 which are 
in the different configurations of the dipole moments along 
[100], as shown in Figs. 14(a) and (b). This degeneracy is 
lifted by t c as follows. In Fig. 17, we compare the inter-layer 
charge fluctuation process by t c in COiy3 with that in COi/ 6 . 
The energy reduction due to this type of the charge fluctuation 
is of the order of t^/^cNNN; this fluctuation easily occurs, 
because the intermediate-state energy Knnn is the smallest 
Coulomb interaction in the model of Hy t - In the CO^ struc- 
ture, the charge fluctuation occurs at the sites represented by 
the black (white) circles connected by arrows in the upper 
(lower) layer in Fig. 17(a). On the other hand, in the case of 
the CO1/6 structure, the number of sites where this type of 
charge fluctuation occurs are half of that in CO1/3, as shown 
in Fig. 17(b). In other CO's, i.e. CO^'s and COi/±, this fluc- 
tuation process is prohibited. 

In order to support this scenario, we calculate the energies 
by the second-order perturbation in terms of i a b and t c . A gen- 
eral form of the energies in the CO's, shown in Fig. 16(a), are 
given by E/V ahNN = 2 - a(t ah /V ahNN ) 2 - /3(t c /V ahNN ) 2 



where a and /3 are positive numerical constants. The sec- 
ond and third terms imply energy gains from the intra-layer 
and inter-layer charge fluctuations, respectively. We have 
(a, (3) =(2.00,2.76) in polar C0 1/3 , (2.60,1.83) in C0 1/6 , 
(2.67,0.67) in C0 1/2 -A and (2.67,1.33) in C0 1/4 . That is, in 
the polar COi/ 3 , a is smallest, while /3 is largest among them. 
This indicates that the stabilization of the polar CO1/3 phase 
is mainly caused by the inter-layer charge fluctuation due to 
t c . 

5. Summary and Discussion 

In this section, we discuss the relation between the present 
results and our previous results in Ref., 27 - 1 and their im- 
plications for the layered iron oxides. In our previous ex- 
aminations, we study the dielectric properties in a paired- 
triangular lattice; the long-range Coulomb interactions are in- 
troduced between the electronic charges which are treated as 
classical variables. The model Hamiltonian was analyzed by 
the mean-field method and the classical Monte-Carlo sim- 
ulation in finite size clusters. It was shown that the three- 
fold charge ordered state, i.e. CO^, does not appear at 
zero temperature except for the fully frustrated point, corre- 
sponding to KNNN/KbNN = 0.6 in Fig. 7 in Ref. 27) and 
VcNNN/VkbNN = 0.5 in Fig. 15, and is stabilized in finite tem- 
peratures due to thermal fluctuation. The obtained V-shape 
phase diagram in the temperature and Vcnnn/Kann plane 
(see Fig. 7 in Ref. 27) ) is similar to that in Fig. 15 where the 
electron transfer integral corresponds to the temperature. This 
comparison suggests that the thermal and quantum fluctua- 
tions play similar roles for stabilization of the polar three-fold 
charge ordered state. Naively speaking, the CO^ state is not 
a robust charge ordered state, since, as shown in Fig. 14(a), 
the intra-layer NN Coulomb interactions at black (white) cir- 
cles in the upper (lower) plane are canceled out. Therefore, 
the thermal motion and the quantum virtual motion of charges 
at these sites easily occur, and these contribute to the entropy 
gain and the transfer energy gain in the classical and quantum 
cases, respectively. This is the so-called order-by-fluctuation 
process. 

One notable difference between the two types of the fluctu- 
ation is that the polarization correlation P defined in Eq. (13) 
is weaken with increasing system size in the finite temper- 
ature classical calculation, while it is robust with increasing 
L in the present quantum calculation at T = 0. In the clas- 
sical Monte Carlo simulation in finite temperatures, a large 
number of non-polar structure coexists with the polar COif 3 
structure. As a result, the thermal expectation value of the po- 
larization correlation tends to decrease with the system size. 
On the other hand, in the present VMC calculation where the 
long-rage charge ordered states are assumed in the variational 
wave function, the fully polarized charge ordered state is en- 
ergetically stabilized due to the inter-layer electron transfer 
and the realistic crystal structure, as explained above. 

Through these two theoretical examinations, i.e. the classi- 
cal Mote-Carlo simulation in finite temperature and the VMC 
simulation at T = 0, we conclude that charge fluctuation 
plays essential roles in the dielectric properties in layered iron 
oxides. These theoretical results have implications for the ex- 
perimental results in layered iron oxides. The charge fluctua- 
tion between Fe 2+ and Fe 3+ has been observed in the optical 
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absorption experiments and the Mossbauer spectroscopy mea- 
surements. 31 ' 4 ^ It was suggested in the temperature depen- 
dence of the optical absorption spectroscopy that the charge 
fluctuation remains to be large even far below the charge or- 
dering temperature. This is consistent with our scenario that 
the thermal/quantum charge fluctuation stabilizes the polar 
charge order. The diffusive nature in the dielectric constant 
observed near the charge ordering temperature is also related 
to the remarkable charge fluctuation predicted in the COi/ a 
state. Another relation of the present results to the layered iron 
oxides is the several different charge ordered states observed 
by changing the rare-earth ion in i?Fe2C>4. In LuFe204, the 
three-fold charge order associated is stabilized below 320K. 
On the other hand, in YFe2C>4, a sequential phase transition 
occurs as (the charge disordered state) — s- (the polar three-fold 
CO) —5- (the non-polar CO's with different periodicities) by 
decreasing temperature. 41 ' That is, the polar CO phase only 
appears in a middle range of temperature. This result is con- 
sistent with our results of Fig. 15 and Fig. 7 in our previous 
paper, 27 ' when we speculate that a change of the R ion cor- 
responds to a change of the Coulomb interaction parameter 
VcNNisr/VabNN and/or the electron transfer intensity. 

In summary, we examine roles of the quantum fluctuation, 
i.e. the electron transfer effect, in the electronic ferroelectric - 
ity. In particular, we focus on the dielectric properties in a 
paired triangular lattice, motivated by the multiferroic lay- 
ered iron oxides. The variational Monte-Carlo simulation is 
applied to the three types of the extended V — t models. It is 
found that the quantum transfer between the triangular lattice 
layers tends to promote the three-fold polar charge ordered 
state. Both the spin degree of freedom and the realistic crystal 
structure in the layered iron oxides reinforce the electric po- 
larization. Present roles of the quantum fluctuation on the fer- 
roelectric transition are in highly contrast to the conventional 
manner of the quantum dielectric fluctuation in the hydrogen- 
bond type ferroelectricity and the quantum paraelectric ox- 
ides. 
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